Spiral phase and phase separation of the double exchange model in the large-S limit 
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The phase diagram of the double exchange model is studied in the large-S limit at zero temperature 
in two and three dimensions. We find that the spiral state has lower energy than the canted 
antiferromagnetic state in the region between the antiferromagnetic phase and the ferromagnetic 
phase. At small doping, the spiral phase is unstable against phase separation due to its negative 
compressibility. When the Hund coupling is small, the system separates into spiral regions and 
antiferromagnetic regions. When the Ifund coupling is large, the spiral phase disappears completely 
and the system separates into ferromagnetic regions and antiferromagnetic regions. 



The colossal magnetoresistance efFect(CMR) was discovered in hole-doped manganese oxides such as Lai-xSr^MnOa 
and Lai^xCaxMnOs P|. Various experiments have revealed that these materials have very rich phase diagrams 0- In 
these materials, the three t2g electrons form a localized S = 3/2 Kondo spin at each manganese site and Cg electrons 
form a conduction band. The degeneracy of the two orbitals is lifted by Jahn- Teller effect and only one of the 
orbitals is close to the fermi energy. The eg spins interact with t2g spins through Hund coupling. The single orbital 
double exchange (DE) modelQ is the simplest description of this system. There are other important factors in this 
system, such as Jahn- Teller phonons, exchange interaction between the Kondo spins, and Columb interaction. It is 
important to study the DE model and identify its role in the CMR systems. 

The ferromagnetic phase of the DE model have been studied in the past 0, 0| . De Gennes studied the DE model 
with the exchange interaction and found that the system goes into a canted antiferromagnetic (CAE) phase near 
half filling ^ . However numerical studies P| discovered phase separation in the simple DE model near half filling in 
one, two and infinite dimensions. There is some evidence that it occurs in three dimensions as well. In addition, an 
incommensurate phase was found when Hund coupling is relatively small. Phase separation was also found in DE 
models with exchange interaction and with Jahn- Teller phonons 0. 

Some theoretical studies 8', ^ ll0lllTl ll2] have been focused on phase separation in the DE model. Phase separation 
was also found in the limit of infinite Hund coupling near Curie temperature, but absent in zero temperature 
The DE model was mapped onto & t ~ J model in Ref. where phase separation was found not only near half 
filling ^3 , but also near zero filling. The stability of canted antiferromagnetic state was examined in Ref. ^3 and 
phase separation was found in certain parameter regions, but it is unclear that without exchange interaction whether 
the phase separation exists or not. In Ref. and Ref. phase separation was found in the DE model within the 
dynamic mean field approximation. 

In this paper, we study the zero-temperature phase diagram of the simple DE model in the large-S limit in cubic 
and square lattices. At zero doping, the system is in an antiferromagnetic phase. At a high level of doping, the 
system goes into a ferromagnetic phase. We find that in the region between the antiferromagnetic phase and the 
ferromagnetic phase, the spiral state has the lower energy than the CAF state. However, the spiral state is always 
unstable near half filling and subject to phase separation. When the ratio of Hund coupling to hopping energy is below 
certain critical value, the separation is between the antiferromagnetic phase and the spiral phase. Above this ratio, 
the spiral phase is always unstable and the system always separates into antiferromagnetic phase and ferromagnetic 
phase. There is one qualitative difference between 2D and 3D systems. In 3D, in the limit of zero Hund coupling, 
the transition between the ferromagnetic phase and the spiral phase occurs at hole density around 0.55; in the same 
limit, the critical density approaches 1 in 2D. 

The Hamiltonian of the DE model is given by 

n = -tY,clc,,~jY,s,-s,, (1) 

where is the electron spin with total spin 1/2 and is the Kondo spin with total spin S. In this paper, we consider 
the Kondo spins as classical spins which is exact in the large-S* limit. The quantum correction to this approximation 
is smaller by a factor of 1/ S and vanishes in the large-5 limit. The classical spin approximation is also equivalent to 
trial wave- functions with static spin configuration. Thus it can provide an upper bound for ground state energy. 
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We consider the possible phases of the system. For convenience, we use the electron operators fi and hi which have 
spins parallel and opposite to the Kondo spin, 

QT = cos(0,/2)e'"'/, + sin(0,/2)e-''("-+^')/i„ 
= -cos(0,/2)e-^"'/i,+sin(0,/2)e^("-+^')/„ 

where ai is an arbitrary phase factor, and 9i and 0^ are the azithmuth and polar angles of S;. For the Kondo spins, it is 
convenient to use the Schwinger-boson representation IQj, 6i| = -^25 cos(6'i/2)e*"' and bii — \/2S' sin(6'i/2)e**^'^'+"'-' 
with the constraint b*^bi^ — 2S. In the classical spin approximation, the Schwinger-boson operators are just 
complex numbers. In terms of the new fermion operators and the Scwinger-boson operators, the Hamiltonian is given 

by 



"—is 



+ h]h,) ^ bib,, + (flh,{b*^^b*^ - b;^b*^ + h.. 



-^E(/'/^-'^l'^')- (2) 

i 

In the ferromagnetic phase, all the Kondo spins are aligned in the same direction, and the electron hopping amplitude 
is diagonal and at maximum, b^bja = 23 and bi^bji — bj^bii = 0. The h- and /-fermions are free particles with the 

dispersion given by i?k = Ck ± \ JS, where ek = — 2t cos fcaO, a is the lattice spacing and d is the dimensionality. 

When the electron bandwidth Adt is bigger than the Hund energy splitting JS, for less than half filling, only the 
lower band is occupied and this fully-magnetized state is an exact eigenstate the Hamiltonian. When JS < Adt 
and sufficiently close to half filling, both bands become occupied and the ferromagnetic state is no longer an exact 
eigenstate but an approximation. It is straight forward to go beyond classical spin approximation and show that the 
spin fluctuations are essentially described spin waves with quadratic dispersion 0,0], much like in a typical Heisenberg 
ferromagnet. At zero temperature, the spin waves are frozen out and do not provide any significant change to the 
classical spin configuration. 

We have also considered the possibility of a totally disordered state. However, in the classical spin approximation 
it is insufficient to consider such a state. We use a mean-field Schwinger-boson formalism '16"] instead. But we found 
that the saddle point of the disordered state does not exist at zero temperature for any positive integer or half-integer 
Kondo spin in dimensions higher than or equal to two. 

Exactly at half filling, the system is in an antiferromagnetic phase. For simplicity, we choose the Kondo spins 
aligning in the x-direction, 6i| — -v/S'e"**^'^' and bii ~ x/S'e'*^'^', where Q = (5^,^,5^)- The hopping matrix 
becomes completely off-diagonal, X^cr ^io-^j'^ ~ ^ ^'^^ ^ ^jT^U ^ —215" sin (Q • (R^ — Rj)). The fermions form 

two bands with dispersion given by i?k = '^^4^ ' + where ej^ — —2t sin k^a. The free-fermion operators are 
linear combinations of the /- and /i-operators. At half filling, only the lower band is occupied. 

In comparison, when JS > Adt, the lower band of the ferromagnetic state are completely filled and the total hopping 
energy is zero at half filling. It is clear that the antiferromagnetic phase has lower energy in this case. When JS < 4dt, 
the total energy difference between these two states is given by 

AE = £'j^iljS)^+e^- ljS)N{e)de + {^JS)^ + - e)Nie)de, (3) 

where N{e) = ^j^. (5(ek — — J2k ^i^w ^ The ferromagnetic state has higher energy in this case as well. 

It is not difficult to go beyond the classical spin approximation to show that there are linearly dispersed spin waves 
in the antiferromagnetic state. Here the quantum fluctuations are stronger because the classical antiferromagnetic 
spin state is not an eigenstate of the Hamiltonian. The actual staggered magnetization is smaller due to fluctuations. 
But it does not change the fact that the antiferromagnetic state has lower energy since the true antiferromagnetic 
state has even lower energy than the classical state. 

The properties of ferromagnetic state are totally different from those of the antiferromagnetic state. It is not a 
total surprise that other states have lower energy in the region between the two phases, as found in Ref. 0]. Here we 
consider two homogeneous states which are natural candidates, the spiral state and the CAF state. 

The spiral state was considered as a possible description of doped high Tc systems 0|, although few supportive 
experimental evidence have been found so far. The spiral state is a compromise between the hopping of the holes and 
the effective superexchange interaction of the electrons. The ferromagnetic and antiferromagnetic phases are just two 
extreme cases of the spiral phase. The possibility of a spiral phase was considered in the DE model with exchange 
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interaction in the large Hund coupling limit [l5l | . It was found that there is a transition from the ferromagnetic phase 
to the spiral phase near zero doping as a function the of exchange coupling. However, it was not addressed whether the 
spiral phase exists or not in the simple DE model with arbitrary Hund coupling and without any exchange interaction. 

In a spiral state, the Kondo spin is aligned in a certain plane and the angle of the spin is a linear function of its 
site position vector. For simplicity, wc consider the spins in the x-y plane and choose the classical Schwinger-boson 
operators to be bi-\ = -v/S'e and bii = y/Se^'^'^' . The Hamiltonian is now given by 



JS 



■Y.^flh-h\h.)- (4) 



For symmetry reasons, we consider the Q-wavevector in (1,1,1) direction which is equivalent to all the other 
(±1,±1,±1) directions. After diagonalizing this Hamiltonian, wc obtain two free fermion bands with dispersion 
given by 



E'k = cos(Q^a)ek ± \l — h sm{Q ^a)'^ e'l- (5) 

The ferromagnetic phase is a special case with Q = and the antiferromagnetic phase is another special case with 
Qx — tt/2. The wavevector Q of the preferred spiral state can be obtained by minimizing the total energy. 

In a CAF state, each sublattice has different magnetization directions. The uniform magnetization coexists with the 
staggered magnetization. The ferromagnetic and antiferromagnetic states are also two special cases of the CAF states. 
To study the CAF state, we take the z-direction to be the direction of the uniform component and cc-direction to be 
the direction of the staggered component, with the ^-values given by 6^]- = bi and bii = (— 1)'62. The Hamiltonian 
can be easily diagnalized. The fermions form two bands with dispersion given by 



E^ = ±^Cl^r + el + JS^e^, (6) 

where the uniform magnetization is given by Sz — — |&2p)/2. 

We compared the energies of different states. The CAF states always have higher energy than the spiral state. 
In Fig. nj the energy of the 3D DE model is plotted as a function of hole density. As shown in this figure, the 
ferromagnetic phase always have the lowest energy when the hole density is above certain value; below this density, 
the spiral phase has lower energy. At half filling, the antiferromagnetic phase has the lowest energy. However, near 
half filling, as shown in Fig.^ the spiral phase is always unstable due to its negative compressibility and it is subject 
to phase separation. One of the phases in phase separation is the antiferromagnetic phase. Depending on the value 
of JS/4dt, the other phase can be either ferromagnetic phase or the spiral phase. 

To determine the phase boundary where phase separation occurs we use in the following method. Let x be the 
overall hole density of the system and y be the ratio of the volume of the ferromagnetic regions or the spiral regions 
to the total system volume, the total energy of the system per site is given by 

Eix,y)^il~y)EA+yEB{-), (7) 

y 

where Ea is the energy of the antiferromagnetic phase per site and Eb{x) is the energy of the ferromagnetic or the 
spiral phase per site at hole density x. Since the energy y) is at minimum as a function of y, the stability 
condition dyE{x,y) — produces the equation 

Ea^Eb{-)--E'b{-). (8) 

y y y 

The solution of eq. ijSJ) is given by 

- = x„ (9) 

y 



where the constant Xc is the critical hole density of the system since the phase boundary is determined hy y = 1. The 
hole density of the spiral regions or ferromagnetic regions x/y is always equal to the critical hole density Xc- 
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We find that the system separates into the spiral phase and the antiferromagnetic phase when ^ is smaller than 
certain value. When ^ is bigger than this critical value, the spiral phase disappears completely and the phase 
separation is between the ferromagnetic phase and the antiferromagnteic phase. In 3D, the critical ratio is about 2/3; 
it is about 0.84 in 2D. In both cases, the critical hole density of phase separation boundary vanishes in the limit of 
J — > oo and also in the limit of J — > 0. 

The complete 3D and 2D phase diagrams are shown in Fig. |21 and Fig. 13 Overall there are five phases: the 
antiferromagnetic phase, the ferromagnetic phase, the spiral phase, phase separation between the antiferromagnetic 
phase and the ferromgnetic phase, and phase separation between the antiferromagnetic phase and the spiral phase. 
The phase transitions are likely to be second order because both the energy derivative and the sizes of various phase 
regions are continuous across the phase transition lines. The 3D phase diagram has a tricritical point at « |, 
X ~ 0.38; the tricritical point of 2D phase diagram is located at ^ « 0.84, x « 0.35. However, the 3D and 2D phase 
diagrams has one major difference. In 3D, in the limit J — )■ 0, the critical density between the ferromagnetic phase 
and the spiral phase is 0.55; in 2D, it is unity. 

In Fig. 131 numerical data from Ref. ;€] are also plotted. The agreement is reasonably good considering that the spin 
fluctuations could quantitatively modify our results and that the finite size effect could affect the numerical results. 
Although it is beyond the scope of our paper, we expect that the spin fluctuations of the spiral states are described 
by quadratically dispersed spin waves, similar to the ferromagnetic case. The dispersion should becomes more linear 
as doping decreases. It is very unlikely that the spin waves will make any qualitative changes to the phase diagrams 
in 3D and 2D, although more accurate numerical phase diagrams are also needed to make comparison. 

In one dimension however, quantum fluctuations in principle destroy any broken symmetry states. But as found 
in Ref. 6], phase diagrams of ID systems are very similar to 2D phase diagrams. It is important to go beyond the 
classical spin approximation to study ID systems and the effects of spin fluctuations. 

In conclusion, we flnd that in the large-S limit of the DE model the spiral phase is an intermediate phase between 
the ferromagnetic phase and the antiferromagnetic phase when the Hund coupling is relatively small. Near half flUing, 
the system subjects to phase separation. The 2D phase diagram that we have found is consistent with the numerical 
results y. Future studies on more sophisticated DE model is needed to relate to CMR systems. 
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FIG. 1: The energies of the 3D DE model in the large-S limit as a function of hole density x for (a) JS = Qt and (b) JS = 19>t. 
E{x) is the total energy per site subtracted by chemical potential — i(l — x)JS. The solid lines are the lowest energies of the 
system. In (a), the dotted line is the energy of the spiral state in the phase separation region. The energy of phase separation 
is given by the straight line below it. The star marks the transition point from the ferromagnetic phase to the spiral phase. In 
(b), the dotted line is the energy of the spiral state in the phase separation region. The star marks the hypothetical transition 
point between the spiral phase and the ferromagnetic phase. 
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FIG. 2: The zero temperature phase diagram of double exchange model in cubic lattice in the large-S limit. The ferromagnetic 
phase (FM) exists at high hole density and the antiferromagnetic phase ( AFM) exists at half filling. When < | , the spiral 
phase (SP) appears at low hole density. Phase separation between SP and AFM occurs near zero doping. When > | , phase 
separation between FM and AFM appears at low doping. 
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FIG. 3: The zero temperature phase diagram of double exchange model in square lattice in the large-S limit. The black spots 
are marks of of the phase boundaries taken from Ref . [6] . 



